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Abstract 

The SU(N) generalization of the multi-channel Kondo model with arbitrary rectangu- 
lar impurity representations is considered by means of the Bethe Ansatz. The thermody- 
namics of the model is analyzed by introducing modified fusion equations for the impurity, 
leading to a simple description of the different IR fixed points of the theory. The entropy 
at zero temperature is discussed; in particular the overscreened case is explained in terms 
of quantum group representation. 
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1. Introduction and basic equations 

The Kondo model with antiferromagnetic coupling was studied in 1964 [1] and found 
account for the low-temperature resistivity of metals in which magnetic impurities were 
present (Kondo effect). The perturbative calculations led to divergences as the temperature 
was lowered, indicating an inaccessible strong-coupling IR region. This stimulated great 
theoretical interest in the field. Several approaches were used to study the Kondo effect in 
more detail ; here we are particularly interested in the fact that the model, formulated 
in terms of one dimensional fields, was found to be solvable by Bethe Ansatz techniques 

BE)- 

In this paper we study the solution of a generalized Kondo model, describing fermions 
carrying color (spin) and flavor (channel) in the fundamental representations of SU(N) 
and SU(f) respectively, interacting with a localized impurity. The impurity is in a (n x 
/) rectangular representation of the color, and in a trivial representation of the flavor 
group. The model combines and generalizes the multichannel Kondo Model of Nozieres 
and Blandin [[7] applicable when the impurity orbital structure is taken into account, and 
the SU(N) version of Coqblin and Schrieffer || valid for rare earth materials with strong 
j — j coupling. 

We shall diagonalize the model (section 1); the Bethe-Ansatz basis separates in a 
natural way into charge, spin, and flavor sectors. This separation, when the interaction 
with the impurity is turned off, corresponds to the decomposition of the free fermion CFT 
describing the electrons into a U(l) sector, a SU(N) level / WZW and a SU(f) level 
N WZW. Only the spin sector interacts non-trivially with the impurity as we shall see 
by studying the associated Bethe Ansatz Equations. In this sector, we shall find that 
two symmetries play a key role: the original SU(N) symmetry of the whole theory, and 
a new hidden quantum group SU q (N) symmetry where q = exp(2z7r/(/ + N)). These 
symmetries will naturally appear in the structure of the effective Bethe Ansatz Equations 
describing this sector. It should be noted that we are here extending the well-known 
connection between the WZW SU(N) level / theory, i.e. the non-interacting spin sector 
of the electrons, with quantum group SU q (N) - through the braiding matrices or the 
truncated fusion rules - to our interacting (non conformal) theory: the fact that the 
quantum group symmetry is preserved is probably equivalent to the fact that integrability 
is preserved. 

We shall next briefly study the ground state and the low-lying excitations (section 
2), and then move on to the thermodynamics of the system. In order to do so we shall 
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consider a scaling limit procedure where the linear density of the electrons D, which can 
be thought of the order of the depth of the Fermi sea, plays the role of a UV cutoff of 
the theory and is sent to infinity while keeping the physical energy scale, De~ 27T / Nc , finite. 
We can then introduce fusion equations dressed by finite temperature for the impurity; we 
quote equations (3.19): 

X?(C + ™/N) X ){C - in/N) = xT\OxT\0 + ^ +1 (C)x, r -i(C)e- 2 ^ sin ^ N ^ . (1.1) 

Here ( is a parametrization of the temperature (T oc e _< >) and x r j is the contribution 
of the impurity whose representation is a rectangular rxj Young tableau to the (finite 
temperature) partition function. So this system of equations connects impurities with 
differents spins (intuitively impurities with high spin are made out of smaller spin impurities 
through the fusion procedure), and it contains all the physics of the model. Indeed we shall 
use it to compute physically interesting quantities in the low and high temperature regimes 
(section 4) . Finally we shall tackle with the interpretation of the zero temperature entropy 
(section 5) using quantum group arguments. 

1.1. The model. 

The Hamiltonian to be studied is, 

n = -ij $ttx)d x i>ttx)dx + JC(o)<V\ 6 (o)r^ V- (1-2) 

The energy spectrum of the electrons is taken to be linear in the momentum as we 
shall be studying universal properties near the Fermi surface. 

The electrons are in the fundamental representation a a (corresponding to Young 
tableau consisting of one box), and the impurity in the representation ta, given by a 
rectangular Young tableau with / columns and n rows, 1 < n < N — 1. The impurity is a 
singlet of the flavor group SU(f), while the electrons are in its fundamental representation. 
The fermionic field ipf(x) annihilates an electron at x with spin (or color) index a, a = 
1, N and flavor index i,i = l, Note that from the point of view of 2D quantum field 
theory, there is only one chirality of electrons ("right-movers"). The operator X q ' t 'a^X^ 
represents the impurity spin operator in a representation {r^}, where the impurity field 
X a is taken to be fermionic and subject to the constraint ^ Q X a X a = 1- Summation over 
all indices is implied, A = 1 . . . iV 2 — 1 being a 5u(iV) Lie algebra index. 
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The interaction with the impurity breaks the U(fN) symmetry of the free hamiltonian 
down to U(l) char9e x SU(N) spin x SU(f) flavor . It will be implicitly assumed that all the 
flavor levels are equally populated. 

We shall find that the model possesses a variety of IR fixed points, whose nature 
depends on the symmetry structure in the flavor sector and on the spin representation 
n x I, generalizing the familiar N = 2 case (the multichannel Kondo model P-|TT|). We 
shall identify the mechanism underlying the appearance of these fixed points as dynamical 
fusion by which the electrons form spin complexes whose interaction with the impurity 
leads to a new behavior in the infrared |T(J. These complexes consist of / electrons fused 
into local objects that transform according to Young Tableaux of one row of length /. 
These composites interact with the impurity and determine low energy properties of the 
model. 

In the Bethe-Ansatz approach a precise description of the formation of these compos- 
ites can be given. To do so, a careful cut-off procedure needs to be introduced to allow 
the formation of the composites while maintaining integrability in the presence of a finite 
cut-off. The linearized hamiltonian propagates separately the charge-spin-flavor degrees of 
freedom that make up the electron. Therefore the effect of flavor on the spin degrees of 
freedom is is recovered only in the full space. To follow the dynamic coupling of spin and 
flavor we add some curvature which maintains the identity of the electron while allowing 
its components to interact. In the end of the calculation the cut-off is sent to infinity 
and the curvature removed. Already in the free field theory the resulting theory is quite 
involved, and even the counting of states is not trivial [|12]. However, the resulting basis 
is the natural one in which to turn on the impurity, it is the zero order approximation in 
the sense of degenerate perturbation theory. 

The scheme consists of the following elements: 

o A second derivative term with a cutoff, A, 



This term explicitly provides an energy cut-off. Furthermore, it introduces curvature 
into the electronic spectrum and breaks charge-spin-flavor (CSF) separation. Once 
the electron composites are formed, and the low-energy spectrum of the theory is 
identified, the cutoff is taken to infinity. 

Adding the term ( |1.3| ) also imposes restrictions on the form of the eigenstates which 




(1.3) 
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can be expressed in terms of the following counterterms without which the model is 
not integrable for finite A. 
❖ An electron-electron interaction term, of the form 



2J J ^(x)^(x^(x)^(x)dx, (1.4) 

The term has no effect on the spectrum once the cut off is removed and no impurity 
is present, independently of the value of the coupling J. The linearized spectrum has 
a large degeneracy which is removed when the interaction with the impurity is added. 



The addition of ( |1.3|) and (|1.4|) provides a way of finding the eigenstates, as we will 
show below, 
o A counterterm 7i cc , of the form 

Hcc = il rt( x ) v ( x )^( x ) dx > 

with 

V(x) = ||t (S'(x +0 ) + 8'(x-°)), (1.6) 

needs to be added to the Hamiltonian in order to preserve integrability at the origin; 
this term vanishes once the cutoff is removed, and plays no further role in the problem. 
Eigenstates of ( |1.2|) with N e electrons and one impurity are of the form 

\f>= J x\<h, /••/;; ( ^({.r^)\;,n '•;?(•'•,)<' >• a-?) 

3 3=1 

where the fermionic field Xb creates the impurity at x = 0. The amplitude F satifies the 
differential equation h\F >= E\F >, where the first quantized Hamiltonian h takes the 
form 

N 



h = ^{-idj - + 2JS(x J )a A r A } 



2A 

3 1 (1-8) 

with Pji(Vji) being the spin (flavor) exchange operator. When the impurity is also in the 
fundamental representation, we can write 

h = Y,{-idj - (A _1 )9| + 2J8( Xj )P j0 ) 

3=1 (1-9) 



Kj 3=1 
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We see that the interaction terms act only when electrons coincide at the same point 
or at the impurity site. Hence, the eigenstate amplitudes are combinations of plane waves 
with pseudo-momenta kj , and have coefficients that depend on the ordering of the electrons 
and the spin and flavor indices. These coefficients are related through products of electron- 
electron and electron-impurity S'-matrices that we will write below. Here we will only write 
explicitely the results for an impurity in the fundamental representation. 

The electron-impurity S'-matrix can be written, to first order in 1/A, 

\ Xj + l-ic J 

where 

x i' 1 + j2 \ k o 2J 



In the scaling limit, J and c have the same scaling behavior. Notice that ( |1.10| ) is trivial 
in the flavor sector. 

The electron-electron S'-matrix is of the form 

Xj - Xi - icPji Xj -Xi + icVji 

bjl = — r : : 7 7 — —. — ■ \ lA2 -) 

Xj — Xi — ic Xj — Xi + ic 

if we set 

J =IT^ (L13) 

Integrability is guaranteed since Sjo and Sji satisfy the Yang-Baxter conditions 

Q Q Q C Q Q 

JjlJjOJIO — JlO^jO^jl 

SjiSjkSik = SikSjkSji 
Finally, the energy eigenvalue of a N e electron state is of the form 

^ = 5>(1 + ^)- (1-15) 

3=1 



(1.14) 



In order to determine the spectrum, we impose periodic boundary conditions, and solve 
the corresponding eigenvalue problem. The procedure is standard [|13| and we skip here the 
details. The result is contained in the Bethe Ansatz Equations (B.A.E.) which we proceed 
to write down. Each of the degrees of freedom - charge, spin and flavor - is described by 
a set of variables whose number depends on the symmetry of the particular state. The 
charge degrees of freedom are given by the set {kj, j = l,...,A^ e }. The spin degrees of 
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freedom are parametrized by the sets {A 7 , 7=1,..., M r , r = 1, . . . , N — 1}. Finally, the 
flavor degrees of freedom are represented by the sets {u; 7 , 7 = 1, . . . , M r , r = 1, ...,/ — 1}. 
The set of integers {M r , r=l,...,iV — 1} specify the symmetry of the spin component 
of the wave function. Similarly, the quantum numbers {M r } specify the symmetry of the 
flavor component. 

The equations are: (L size of the periodic space) 

kjL = Ai-q + A.Q+zf fC ^-(l + A.Q-zf 
llAi-(l + A,)-z|ll^-(l + A,) + zf' 

M 1 ,1 , ,1 , • N e 1 . , n . . c M 2 , ,1 , ,2 



TT ~7 ^ 1 _ = TT ^7 ^ 1T/, ^ T '2 TT 
11 a/} - o/J - ic na;i-(l + A,)-z| llo;i-^ -if 

n : ? = n n : ? - 2 ,...,/-i, 

£ ^-Ai+»c _fi A^-(l + A 3 ) T if 5! A 7 -A| + zf 
■llAi-Aj-fc llAi-(l + A,)-z§llAi-A|-zf' 



-7 


-<4 


— zc 






+ ic 


tor 




— ic 


A 7 




+ ic 


A 7 




— ic 


A n 

7 


— A n 


+ ic 


A™ 


-Aj§ 


— ic 


A r 

7 




+ ic 


A r 

7 


— A r 


— ic 



(1.16) 



M n An A n 1 „•„ A n i Al c M \n A t , a c 

A 7 -Mfj TT TT A 7 y V +Z 2 
■ - . - ■ \n_A]C 11 11 A n _ At _ -c ' 

_„ n n A ?: A f! i >-= 2 w-w». 

0=1 to t=r±1/3=1 A 7 A /3 Z 2 

The next step is to solve the equations for all possible states, identify the ground state 
and the low energy excitations above it. Subsequently, by summing over all excitation 
energies we shall obtain the partition function. 

The B.A.E. are a function of the cutoff A which eventually is sent to infinity. In this 
limit the equations reduce to a smaller set once the correct ground state has been identified. 
It corresponds to string solutions (see below) leading to electron composites which interact 
most efficiently with the impurity. The ground state and low lying excitations lie in a 
sector of the theory given by solutions of a particular form - f -strings. Solutions of this 
type are SU(f) flavor singlets which allows them to have maximally large SU(N) spin. 
We shall see that this class of excitations is characterized by a scale T = De~^. When 
strings are broken to form flavored excitations we expect them to be characterized by other 
scales which will tend to infinity as the cut-off is removed and thus do not contribute to 
the impurity dynamics. 
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The formation of composites in flavor corresponds to solutions of the B.A.E. where 
the charge parameters, are complex numbers centered around {^}. Likewise, rank 

r flavor parameters are themselves centered around rank r + 1 solutions |10| . The form of 
the charge parameters is, 

{Xj,j = l... N e } = |p 5 /A + ic U±±- - q) ; q = 1, 2, . . . , f, p 5 real, 5 = 1, iV e //| . 

(1.17) 

while the flavor parameters, 

{u>;, 7 = 1,2, ...,M r }= \p s /A + ij( f ~ r + 1 -q) ; </ = 1, 2, r, 5 = 1, . . . , N e / f 



(1.18) 

where r = 1, ...,/ — 1. These configurations satisfy the B.A.E. in a trivial manner and 
induce fusion in the B.A.E. equations as well as in the form of the wavefunctions. A string 
built on momentum p as its real part induces in the wave function a composite of the form 
exp{ — TjAJ J2j i \ x j ~ x i\ + ip(%i + ••• + x f)} x [•••]) which becomes local as A — > oo. These 
composites will be described by effective fused B.A. Equations obtained by inserting the 
string configurations into the full B.A.E. After removing the cutoff they become 



^ A--A^ + zc *! A; - A^ - ic/2 



;i.i9) 




l + z/c/2\ /A™ + ilc/2 



ifc/2 \A™-ilc/2 



r ~i i 



for each root AL A'e is now the number of composites of / electrons (i.e. originally N e / f). 
( |1.19|) will be our starting point for all subsequent calculations. 
The energy e of a composite of / electrons is given by: 



e iL£ =ii( tt : • ;.v: " i (1.20) 




The (irreducible) representation R of the resulting state is given by the numbers M r 
of roots of type r; indeed, each root A^, moves one box of the Young tableau down from 
row r to row r + 1. In other words one starts with the empty configuration M r = which 
corresponds to the highest possible weight (fig. 1). 



Fig. 1: The highest possible weight of the model. 
We call this representation R . Then 



N-l 



R = R -J2 M 



(1.21) 



r=l 



where a r = e r — e r+1 is the r th simple root of the Dynkin diagram A^-i (e a , a = 1 . . . N 
corresponds to one box in row a of the corresponding Young diagram, see appendix A for 
more details). 

1.2. Continuous B.A.E. 

As we are mainly interested in the thermodynamics of this model, we shall immediately 
write equations in the thermodynamic limit (N e — > oo, keeping the density of electrons 
per unit length D = N e /L and the densities of roots M r /L fixed i.e. of the order of some 
physical energy scale). In this limit, standard calculations lead to continuous B.A.E. Let us 
briefly derive these in a rather formal way which will minimize the amount of calculations. 

We recall that according to the "string hypothesis" , the A^ group into strings of roots 
with the same real part, and a fixed distance of ic between two consecutive roots. 

We therefore introduce some notations. We shall need the function 



0(A) = 1 log A 



ic/2 n f2A\ 

j- = 2 arctan — + tv 

A + zc/2 V c ' 



(1.22) 



where we have specified the determination of the log on the real axis, so that ©(— oo) = 0. 
We also define its "descendants" generated by strings: 



e fel ,..., fcp (A)= e(A + zc(j 1 + --- + Jp )), 



(1.23) 



each j a being an integer (or a half-integer, depending on the parity of k a ) which ranges 
from -(fc -l)/2 to (fc - l)/2. In particular 6 fc (A) = (1/i) log((A - ikc/2)/(A + ike/ 2)). 
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Again, on the real axis, the analytic structure is chosen so that no cut crosses the real axis 
and 0fe 1 ,...,fe J ,(-oo) = 0. 

In the same way we define 

K(A) = — 4-8(A) = — AO ° - - (1.24) 
v ' 2ndA K ' 2ttA 2 + c 2 /4 v ; 

and its descendants if^,...^ . 

If A^. 7 is the center of the 7 th string (1 < 7 < MJ) of length j (j positive integer) 
and of type r {1 <r < N — 1), then by multiplying eq. (|1.19|) for all the j roots of a given 
string and taking the log, we find: 

27r/J ;7 = r)' 1 .V, (-)/.,( A' /: . - 1) + S^QuiA^) 

00 Mfc 00 M* 25) 

-EE - + E E E e ^( A ^ - A %P) ' 

k=l f3=l t=r±lk=l/3=l 

The Ij are half-integers of a given parity. Differentiating once with respect to A and 
introducing the densities of roots and holes and p£ we get: 

~p r j +p r j =6 rl N e K f j(A-l) + 6 rn K lJ (A)-J2K2 !j ,k*P r j + E J2 K ^*Pk (1-26) 

fc t=r±l fc 

where ★ denotes convolution in A space. From now on, unless explicitly stated otherwise, 
all subscript indices are string indices, running from 1 to 00, and all superscript indices 
are Dynkin diagram indices, running from 1 to N — 1. 

p r j is normalized by the condition • j J ^A p'j = J? • i -WJ = • This is slightly 
incorrect since the M r diverge in the thermodynamic limit, only M r /L is physical; but it 
simplifies notations. 



To simplify even further the form of Eq. (|1.26| ) let us introduce the "spectral 
parameter-dependent" Cartan matrices C qr ', Cjk of Aj^-i, and their inverses G qr , 

For 1 < r < N — 1, C qr is defined by 

C 9r (A) = <^ r 5(A) - (5 9r+1 + 5 gr - 1 )s(A) (1.27) 

where s is the simple function s(A) = l/(2c cosh(7rA/c)). The same formula holds for Cjk, 
only the boundary conditions 1 < j, k < 00 being different. 
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In order to express G qr and Gjk we use Fourier transform defined by 

4>{k) = J 4>(A) exp(i/c2A/c)dA 

for any function <p(A). We have s(«) = 1/(2 cosh(«)), and 

G^(K) = CT*(k) = 2coth(.) Sinh(( ^7 / ) : r ) f nh(r " ) g > r 

sinh(A^) " ( L2 8) 

Gjk(n) = Gkj(fi) = 2 coth(ft) exp(— sinh(fcft) j > k 

Eq. ( |1.26| ) can be rewritten in a simple manner using these kernels. Let us prove for 
example that K 2 j } k = Gjk for j > k. One starts with Kj(n) = exp(— then adding 
the indices k (resp. 2) amounts to multiplying (for convolution) by smh(kK) / sinh(K) (resp. 
2cos1i(k)), which produces Gjk- 

Thus, Eq. (|1.26|) becomes: 

Pj + Y, C9r * G i'<*f>k = fj (1-29) 

q,k 

where by definition /J = 8 rl N e K fJ (A - 1) + 8 rn K hj {A). Multiplying by C jk we finally 
find: 

Yl C ^ * Pk + ^ * Pi = 5 i /^ rl ^(A - 1) + Sj^siA) (1.30) 

k q 

2. Ground state and low-lying excitations 

In this section, we discuss the nature of the ground state and of the physical excitations 
above the ground state. Before proceeding, the following remark should be made: 

In order to get a simple consistent picture (that sheds some light on the thermody- 
namic results we find in the next sections), we shall intentionally choose to ignore some 
purely discrete effects which are connected with the exact number of electrons. As an ex- 
ample, let us remember that even in the simplest Bethe Ansatz-solvable model, the 577(2) 
spin 1/2 XXX model, the ground state representation depends on the parity of the number 
of spins; but the thermodynamic properties of the model do not depend on it, and one can 
make the simplifying assumption that it is even. In our case, we shall ignore these effects 
by allowing the M r to be sometimes non-integer. Such a procedure should certainly not 
affect the thermodynamic quantities such as the free energy in the L — » oo limit0. 



1 Note that these effects will modify the l/L corrections to the free energy, which are themselves 
of the same order as the free energy of the impurity (cf Eq. (3.3) of next section); but the latter 
should be unaffected at leading order. 
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2.1. Energy. 

According to ( |1.2U| ), the energy (without any magnetic field) is given by: 

M 1 



7=1 



M] 



= -^EE /»i(47- 1 ) 

j 7=1 

= -Dj2 1 dAp}(A)e Lj (A-l) 

j 

(D = N e /L). For future use we shall rewrite this 

where c/J(A) = -D5 rl S f ,j(A - 1). 

The energy can also be expressed in terms of p: 



E = -DJ2J dAC jk *G rl * {-pi + fi)B flj (A - 1) 

r,j,k 

= E g , s ,+Dj2 [ dAp r f (A)G rl * [2 arctan(e 7r ( A - 1) / c ) 

r J 

which we can again rewrite under the form 



E = E g . s . + Y, J dA#(A)#(A) 

r,j 



where explicitly: 
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D5 



if 



2 arctan I tan 



7r N — r 
2 N 



tanh ( -^-(A 



1) + ^ 



N -r 



N 



(2.1) 



(2.2) 



(2.3) 



(2.4) 



(2.5) 



We first give a brief explanation of Q2.3|)-( |2~^|) . E s _ s _ is the energy of the ground state. 
It is obtained for p r j = 0, that is the ground state is filled with the maximum amount 
of /-strings; we shall build it explicitly in 2.3. Above the ground state, excitations are 
created by holes in the sea of /-strings; the other strings do not contribute to the energy. 
Their role should become transparent in next subsection. 
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Let us also see what happens in the scaling limit, when the density of electrons per 
unit length D 1 which will eventually become the UV cutoff of our theory, is sent to infinity; 
then, a typical state of our system will have densities per unit length of /-strings MJ/L 
that will be of order D and will diverge. On the other hand, densities of holes of /-strings 
Mf/L, and of j-strings MJ / L (J ^ /) will remain of the order of the physical energy scale. 
We have already rewritten the energy E — E^ a _ in terms of the only, whose densities 
remain finite. It is interesting do the same thing for the quantum numbers characterizing 
the SU (N) representation of the state. 

2.2. Representation of a state. 

We have seen that a Bethe Ansatz state is in the representation R = Rq — ^ r M r a r , 
where the M r are the numbers of roots of type r: M r = ^2jjMJ. A better way to 
describe the corresponding SU(N) Young tableau is to count the number of columns n r of 
a given length r (1 < r < N — 1). One can easily show the following relation (cf appendix 
A for more details): 

n r = 5 rl NJ + 5 rn l -2j2c qr J M j (2-6) 

Here C qr stands for C qr (K = 0), that is (up to a factor of 2) the usual Cartan matrix of 
SU(N). 

We shall now relate the quantum numbers n r to the numbers of holes MJ = J P r f an d 
of strings Mj (j ^ /). 

We start from equation ( |1.29| ) with j = /, which we integrate from — oo to +oo (i.e. 
take k = 0): 

M r f +J2C qr G ff Mj+ C qr G kf Ml = G ff 5^ + G fl 5- 1 - (2.7) 

where the argument k = is implied for all kernels. Using this formula, one can easily 
compute the sum of eq. (|2.6|). We give the final result: 



'M r f - ^(/-/)2C^Mj ? l<f 



n r = < 



J>f ' q (2.8) 

- E v - ° qrM J + - w rn 1 ^ f 

3>f,q 

The first observation is that we must treat separately the underscreened case (/ > /) and 
the overscreened case (/</): this is directly related to the representation/degeneracy of 
the ground state, which we shall discuss in next subsection. 
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Notice next that the role of the j-strings, j > f, which did not contribute to the 
energy, is now clear: they allow to lower the spin of the system at fixed energy, i.e. fixed 
physical excitations@. This is the usual way, in Bethe Ansatz systems, to select between the 
different irreducible subrepresentations inside the tensor product of the representations of 
the physical excitations (without any j-strings, j > f, the system is in the highest possible 
weight representation). 

On the other hand the j-strings (j < f) do not change the representation of the 
state: in fact they play a role which is quite similar to the j-strings (j > /), but for an 
other quantum number which can be associated with a state, its SU q (N) representation 
(q = exp(2z7r/(/ + N)), so that the set of representations is "restricted to level /"). We 
shall elaborate on this new quantum number in section 5. 

2.3. The ground state. 

According to Eq. (|2.3|) , the ground state is obtained for p r ^ = 0. This means that the 
ground state is filled with a continuous density of /-strings of all types r = 1 . . . N — 1. It 
can then be shown that there are no j-strings (j ^ /); this allows to calculate the densities 



We shall now discuss separately the different cases: 

• / < I (underscreening). We first compute the ground state representation. Applying 
( |2.8|) to the ground state, we find the following Young tableau (fig. 2). 

Let us try to interpret intuitively this Young tableau. First there are a large number 
(of order N e ) of electrons that are not directly interacting with the impurity: they are 
represented by the "trivial" part (i.e. iV rows) of the Young tableau, to the left of the 
impurity. Then there are exactly N — n fused electrons which can be thought as on the 
same site as or glued to the impurity: Their effect is to reduce its spin from n x I to 
n x (I — f) . To confirm this analysis, one can calculate the density of holes of j-strings p^; 



It is known that for a small number of excitations, the complex roots (non-/-strings) do 
not necessarily form exact j-strings: the spacing of their imaginary parts might deviate from the 
string behavior. However this does not modify qualitative analysis. 



p r f from ( |r!Zg| ): 



p r f (A) = N e G rl * S (A - 1) + G rn * (Gff)- 1 * Kif(A) 



(2.9) 



one finds: 



#(A) = 5™(K ltj (A) - G jf * (Gff)~ 1 Kfj(A)) 



(2.10) 
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1 



N 



Fig. 2: Underscreening of the impurity by the electrons. On this example, 
3 fused / = 2 electrons screen the 5x4 impurity, reducing its spin to 5 x 2. 

For j > f, p™ ytz o so there are holes of j-strings; but one can check that there are "not 
enough" to actually create j-strings (this confirms the fact that the ground state cannot 
contain j-strings). The interpretation of this result is the following: the underscreened 
impurity is in a given irreducible representation of SU(N), so its spin cannot be changed. 
However, if one adds one physical excitation into the system, it will create more holes 
of j-strings and allow creation of j -strings: that is so because the tensor product of the 
(non trivial) representation of the impurity and of another non trivial representation is 
necessarily reducible, and the choice between the different irreducible subrepresentations 
is made, as has already been mentioned, by inserting appropriate j-strings (cf Eq. (|2.8|)). 

On the other hand, for j < f, one finds from ( |2.10| ) that p r - = 0: the screened impurity 
possesses a trivial SU q (N) quantum number. 

9 f > I (overscreening) . This time Eq. (|2.8|) implies that the ground state is in the 
trivial representation, which is not what one might naively expect; indeed, by minimizing 
only the interaction part jijjf(0)(T^^ (0)x Q r^x /3 °f the Hamiltonian, one would obtain 
a different ground state representation [^4[ . Let us mention now that what this proves is 
that the RG fixed point J = oo is unstable to the kinetic part of the Hamiltonian and that 
we have a non-trivial IR fixed point J = J*. This will be explained more carefully when 
we study the thermodynamics of the model. 

One can again calculate the densities of holes of j-strings: the formula ( |2.10| ) is still 
valid; but the analysis is reversed. For j > /, one find p^ = 0, which confirms the fact that 
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the ground state is in the trivial representation. On the contrary we now have holes of 
j-strings for j < f: the impurity possesses a non-trivial SU q (N) quantum number, though 
it originally possessed none. Of course its origin is related to the presence of the electrons 
that screen the impurity. We shall discuss this in detail when we consider the entropy at 
zero temperature. 

2.4- Physical excitations. 

We have already mentioned that physical spin excitations are created by inserting 
holes of /-strings. As there are N — 1 types of holes, we conclude that there are N — 1 
types of physical particles labeled by r = 1 . . . N — 1 . The rapidity A at which the hole is 
inserted determines the energy and momentum of the corresponding excitation. We have 
already given in 2.1 the energy e r (A) = g r AA): (Eq. (|2H 



e r (A) = D 



. . TvN-r\ , / 7T a\ N 

2arctan I tan I — — — — I tanh I — — (A — 1) I I + tv- 



2 N \Nc K ) j N 



(2.11) 



As e r (— oo) = 0, these are massless excitations. From general arguments it is clear that 
they have a linear dispersion relation, so that their momentum p r = e r . This can also be 
easily extracted from the Bethe Ansatz equations by writing a phase shift condition on a 
compactified space for one single physical excitation. 

We may at this point go to the scaling limit D — > oo, c —>■ keeping the physical scale 
To = De~ 2ir / Nc fixed. The energy /momentum then takes the simple form: 

c r = „r = 0Tn ■ C*r\ 2nA/Nc 



e r =p r = 2T sm{—)e^" c (2.12) 

which is characteristic of a relativistic (massless) right-moving particle, with 2tt A/ 'Nc its 
rapidity. 

We also know their representation from 2.2: according to Eq. (|2.8[) , inserting a hole of 
type r creates a column of size r in the corresponding Young tableau; that is, the particle 
of type r belongs to the fundamental (or totally antisymmetric) representation of SU(N) 
with r boxes. It can also be shown that they belong to the same representation of the 
quantum group SU q (N). 

These particles interact with each other and with the impurity, which leads to the 
concept of phaseshift. The latter can be extracted from the Bethe Ansatz equations; 
however, we postpone their detailed study to a future publication. 
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3. Thermodynamic Bethe Ansatz equations (T.B.A.) 

We now want to study the model at finite temperature T (T ^> 1/L). In the standard 
way, one derives the T.B.A. by minimizing the free energy F = E — TS with respect to 
p r -\ using ( |1.29| ) one finds that 8S/8p r - equals: {jf- = p^/p^) 

V 

mi + V-) -J2 cqr * * io g (i + (vir 1 ) = | (3.i) 

q,k 



where we have used grj = 6E/8pj, the energy of a j-string of type r (Eq. (|2.2|) ). Note that 
the T.B.A. do not depend on the representation of the impurity. In fact one can say that 
the T.B.A. only describe the electrons (in a way appropriate for studying their interaction 
with the impurity) and not the impurity itself. On the contrary, the fusion equations that 
will be written later describe specifically the impurity. 

3.1. Free energy. 

General T.B.A. formulae imply that 

F = -Tj2j dA/JCAJlogCl + C^CA))" 1 ) 

j,T (3.2) 
(using T.B.A.) =- J2 jdkG* r * C jk * /«(A) [-#(A) + Tlog(l + r£(A))] 

The first term is the ground state energy. In the second term, the explicit expression of 
/| leads to 

F = E'g.g. + LT + FP, (3.3) 
where T is the free energy per unit length of the electrons: 



F = -DTJ2 [ dAG rl (A-l)log(l + ?^(A)) 

r J 

and Fp is the free energy of the impurity: 

FP = -TJ2 J dAG™(A)log(l + i7r(A)). 



(3.4) 



(3.5) 



Here we have introduced a new notation which will prove convenient: G qr (A) = G qr *s(A)] 
G qr has the following Fourier transform: 

= &*(k) = Slnh(( ^ " \ q)K lf^[ rK) q >r (3.6) 
w w sinh(«) smh(ATK) v ' 
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3.2. Scaling limit. 

We now take the limit c — > 0, D — ► oo, keeping To = De~ 2l7 ^ Nc fixed; To is the 
dynamically generated energy scale due to the presence of the impurity. We rewrite the 
T.B.A. equations (371) under the form 



mi + (v-r 1 ) -J2 G9r * c » * mi +vi) = % ( 3 - 7 ) 

q,k 

which become in the scaling limit (Eq. fl2.12p ) 

log(l + (v-r 1 ) -J2 G<lr * C i* * l0 S(! + 4) = 26 jf sin (77) eC (3-8) 

q,k 

where we have successively rescaledi and shifted the rapidity: 

C = ^ + log(T /T). (3.9) 

The shift removes any dependence of the T.B.A. on T. The qualitative picture of 
the behavior of these equations is then the following: there are two asymptotic regimes 
C — > ±00 characterized by limiting values ?7j(±oo). We shall derive these values in the 
next subsection, once we have explained their group-theoretic meaning. 

As the crossover between the two regimes occurs for ( ps that is j^A ~ log(T/To), 
C — > — 00 will be called the high-temperature regime, whereas £ — > +00 will be the low- 
temperature regime. 

Let us also rewrite the free energy in terms of the new variable: we use for example 
G rl (A) = (l/Nc) sin(7rr/A^)/(cosh((27r/iVc)A) - cos(7rr/iV)): 

^=-^E sin (i) S dCe " log(1+?? / ( o) 

r - (3-10) 
FP = -TJ2 / <G™(C-log(T /T)) log(l+i7r(0) 

r J 

The fact that Tq only appears in the free energy of the impurity indicates that it is the 
presence of the impurity which triggered the energy scale dynamical generation. 



3 Note that the rescaling also induces a rescaling of the kernels, according to the rule <f)(Qd£ 
(p(A)dA. 
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3.3. Fusion equations and interpretation of limiting values. 



We shall now work with the T.B.A. equations (|3.8|) , i.e. after the scaling limit has 
been taken. Let us see how these equations are related to fusion equations for rectangular 
Young tableaux. The first step is to analytically continue the functions r] r - to the strip 
| ImCI < 7r /N , so that they possess no zeros or poles for | Im£| < tt/N. Then log(l + 77J) 
and log(l + (Vj) -1 can &l so be continued in the same way. For any function that satisfies 
this analyticity requirement one can then define an inverse operator of G qr which we of 
course denote by C qr ; this operator is not a convolution kernel since it acts in the following 
way: 



C qr * 0(C) = S qr WC + ™/N) + 0(C - in/N)) - (5 qr+1 + s^wc) 



(3.11) 



C qr can be considered as a discrete Laplace operator in (r, () space. As G qr = G qr ★ s, one 
should have C qr = C qr * s: this is indeed the case on condition that one correctly remains 
inside the analyticity strip when performing the convolution integral (fig. 3). 



s(0 








, — r A 


^ : 




- i7t/N 



Fig. 3: s(C) = iV/(47rcosh(iVC/2)) has poles at ( = ±irc/N, so that when 
one considers the convolution kernel s(£ + in/N) + s(£ — in/N), one should 
deform the integration contour so as to remain inside the analyticity strip. 
Then due to 2z7r/A^-antiperiodicity only the pole contribution at C = +i%/N 
remains, so that s(( + in/N) + s(£ — in /N) = 6(C). 

We now take eq. (|3.8|) and act on it with C qr : this annihilates the r.h.s. and we find 
after exponentiation: 



TjUC + iTT/NXiC-iTT/N) 



[l+^+l(0] [l + t/J_!(0] 



(3.12) 



[i + C^CO)- 1 ] [i + lf^C))- 1 ]" 

This system of equations ( "Y-system" ) |15|] is well-known to be related to SU (N) "spectral 



parameter-dependent" fusion equations JT^]. For now, we simply note the "rank- level 
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duality" of these equations: j <-> r, r\ <-> In fact, we have to be careful because of the 
non-trivial r.h.s. of Eq. ( |3.8| ), which does not appear in the framework of normal fusion 
equations. We shall therefore proceed with caution and rederive the fusion equations from 
the very beginning. First we introduce the generalized characters 



X r i (C) = exp 



£G* r *log(l+7/j 



9\ 

L 9 



(3.13) 



(we call them generalized characters because normal characters are the solutions of the 
standard fusion equations, whereas here we deal with modified fusion equations). Notice 
the similarity with Eq. ( |3.5|) or ( |3.10| ): one has 

logx?(0 = -^Fr(T = T e-<) (3.14) 

In fact, as the T.B.A. equations do not depend on the representation of the impurity, one 
can consider that the x r j i f° r an V values of r and j , are related to the free energy FJ of an 
impurity in the representation r x j (r rows, j columns) ! 

One can invert relation ( |3.13| ) using the definition of C qr : 



i^io^^yi?^ (3.i5) 

1 x; +, oxr'c 



(Xj and xf are by convention equal to 1). So far we have not used the T.B.A. (|3.8|) yet. 
We now do so in order to derive "dual" (in the sense of rank-level duality) expressions for 

1 + (Vj)~ X - First 

we introduce a Cjk defined by the same formula ( |3.11| ), and such that 
Cjk * s = Cjk- As G qr ★ Cjk = G qr ★ Cjk, the insertion of the definition of the characters 
in Q3.8p gives: 

log(l + (V-)- 1 ) =J2 d i** l °gXk + 26 jf sin (^) eC ( 3 - 16 ) 

k 

Notice the additional term due to the r.h.s. of ( |3.8| ), which leads to the following modified 
formula: 

1 + (^(C))" 1 = J ? V / V> e 3f Sin{nr/N)e (3.17) 

Xj+i\s)Xj-i\s) 

(Xo = 1) or finally 

r.r(f \ _ Xj+vXjXj-vX) -26 jf sin(7rr7AQe c /o i o\ 

This, with the boundary conditions x°j = xf — 1 an d Xo = 1? allows to compute rjj(() as 
a function of the Xfc(C)- Remembering that x\ ls directly related to the free energy of an 
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impurity in the representation t X k, one can consider that one puts "test" impurities (cor- 
responding to any rectangular Young tableau) in the system in order to find the characters 
Xj, and from there the functions 77J. 

One can now proceed to write down modified fusion equations for the x r j- rewriting 
the trivial identity (|37Lj) = 1 + (|37Lj) we find 

x^(C WAOx-tC-WAO^f^Oxr 1 ^^ ( 3 - 19 ) 



This system of coupled non-linear equations can be represented by fig. 4. 



f-i 



f+i 



1 o 



N-2 



N-l 



Fig. 4: Pictorial representation of the modified fusion equations. The 
marked circles correspond to the position of the modified term. The square 
around the circle n, / means that it is x? which gives the free energy of the 
impurity. 



Just like the T.B.A. equations, the fusion equations contain all the physics of the 
Kondo model. In fact, it is the additional term exp(— 25jf sm(irr / N)e^) which really 
contains all the information on the crossover from the low-temperature regime to the high- 
temperature regime. We shall elaborate on this in next section which contains all the 
calculations of thermodynamic quantities. 

For the moment, in order to see the role of the additional term in a very simple setting, 
let us consider the limits ( — > ±00; then the shifts of ±iir/N of ( become negligible, which 
means that the equations (|3.19| ) reduce to SU (N) tensor product equations for rectangular 
Young tableaux, and the Xj tend to asymptotic values x^(±oo) which should be ordinary 
SU(N) characters. We now compute these values in the absence of magnetic field: 
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• When £ — > — oo, the modification induced by the r.h.s. of ( |3.8|) vanishes and one 
simply obtains ordinary tensor product equations for any j > 1 and r. Therefore 



x r AG) 



(3.20) 



where Xj(G) is the character associated with some matrix G G SU{N) and the rectangular 
Young tableau r X j (cf appendix A for a definition). 

The constraint that 77J > (i.e. x V j > 0) f° r au r an d i, the boundary conditions 
Xo = 1, X-i = an d the asymptotic behavior of as j — > 00 (we shall work out the latter 
in the more general case of the presence of a magnetic field, see 3.4) lead to the unique 
solution G = lsu(N), i- e - the T.B.A. select the Perron-Frobenius eigenvalue of the tensor 
product matrix corresponding to the Young tableau r x j. This gives the following values 

Of T)j(~ °°) : 

Vj(-OO) = Z7Tt — ( 3 - 21 ) 



r(N — r). 

• When Q — > +00, the fusion equations (|3.19|) are cut in two separate parts: for 
j < /, we can consider that we have normal fusion equations with the additional boundary 
condition X r f+i = 0; this, with the usual boundary conditions Xo = 1? X-i = 0? implies 
that 

xJ(+oo) = XjiGf) j<f (3.22) 

where Gf is a SU(N) matrix which can be uniquely determined by imposing 77J > 0: 
(9 = tt/(/ + AT)) 



V 



i(-JV+l)6 



,i(-JV+3)0 



\ 



3 i(AT-3)e 







3 i(AT-l)6> 



(3.23) 



We shall see in section 5 that this character can again be interpreted as the Perron- 
Frobenius eigenvalue of an appropriately truncated (at level /) tensor product matrix. 

For j > f, this time everything happens as if X V f-\ = 0- Together with Xt = 1 (this 
comes out of ( |3.22| )-( j3723|) ) and the asymptotic behavior of x r j as 3 ~^ °°7 ^ reproduces 
the boundary conditions found in the limit £ — ► — 00, except for a shift of / in the string 
index. Therefore: 

X r j(+oo) = x r j-f(lsu(N)) J > f (3.24) 



21 



Finally, 



sinCyfl) sin((j + N)9) 1<j< f 
, sin(r0)sin((iV-r)0) ~ ~ 

r(N — r) 
3-4- Magnetic field. 

In a 5 , L r (A?")-invariant model, the most general magnetic field B one can impose is an 
arbitrary element of the Lie algebra su(N): 

N 2 -l 



B = B aTa (3.26) 



A=l 



where the T4 are the generators of su(iV) acting on the whole Hilbert space of our model. 
Using the SU(N) symmetry one can suppose that B belongs to the Cartan subalgebra; 
this means that B can be described in the fundamental representation of SU(N) as the 
diagonal matrix: 

( Bl ° \ 



B 



(3.27) 



V B N J 

with the condition B 1 +B 2 -\ \-B N = 0. One still has a residual symmetry of permutation 

of the eigenvalues of B which allows to choose: B 1 < B 2 < . . . < B N . This choice will 
prove convenient later. Of course in the SU(2) case B has only two eigenvalues B 1 = —B 
and B 2 = B, and we recover the usual one-component magnetic field. 

The effect of the magnetic field is that, for a given Bethe Ansatz solution belonging to 
an irreducible representation R (characterized by its highest weights, or by the numbers 
M r ), the states with different SU(N) quantum numbers (weights) have a different energy. 
In other words, the corresponding energy level E has been split in several levels; the 
resulting contribution to the partition function then factorizes: Z = ZbZ-j-i, with Zji = 
exp(-E/T) and 

Z B =XR(e- B/T ) (3.28) 

where e~ B l T is the diagonal matrix with eigenvalues e~ fia / T , and xr denotes the character 
of the representation R. Fortunately, this expression simplifies in the thermodynamic limit, 
when one deals with large Young tableaux: then, one can show (see appendix A) that 

_ V-JV-l M r b r /T 

Z B oc e ^r=i 1 (3.29) 



22 



Here one has used the ordering of the B a , and defined b r = B r+1 — B r . Hence, —V 
effectively acts as a chemical potential for the r th type of B.A. roots^L 

The change induced by the addition of a magnetic field in the T.B.A. is straightfor- 



ward: the r.h.s. of (|3.1|) is now: 



f = -^ rl ©/,i(A-l)+J^ (3.30) 

The other equations (|3.2|) - (|3.19| ) are unchanged: the presence of the magnetic field 
is hidden in the asymptotics of 77J as j 1 — > 00. Starting from ( [3.1|) , simple manipulations 
show that 



lim Kj * log(l + n r j+1 ) - K j+1 * log(l + Vj) = % (3-31) 



or for the characters 



lim K 3 * log X r J+ i - K J+1 * log ^ = V ^. (3.32) 

a=l 

4. Calculation of thermodynamic quantities. 

We shall work directly with the Xj an d the fusion equations (|3.19|) they satisfy, since 
they will provide us with the free energy of the impurity, through the correspondence ( |3.14| ) 
that we recall: 

FP(T) = -TlogxHC = log(T /T)). (4.1) 

4-1. High temperature (T 3> Tq). 

The high temperature behavior is governed by the £ — > —00 region of the T.B.A. In 
this region, the modified fusion equations become "free" in the sense that the extra term 
of Eqs. ( |3.19| ) exp(— 25 jf sm(-rrr/N)To/T) = 1 up to 1/T corrections. Let us see what this 
means physically for the impurity: 



When a magnetic field is present, one can show that the solution ( |3.20| ) of the fusion 
equations is still valid, except that the matrix G is no more restricted to be in SU(N) 
(only in SL(N)); in fact, the choice 

G = e- B ' T (4.2) 



Actually, even for B = 0, Zb gives an additional entropic factor in the partition function 
that we have not taken into account so far. But its effect is negligible in the thermodynamic limit: 



only the leading behavior embodied in Eq. ( 3.29| ) matters. 



23 



gives the correct asymptotics (|3.32|) , as can be checked by applying the asymptotic form of 



characters for large Young tableaux (appendix A) to x^, j — ► oo. There are l/£ corrections! 



to the C — > — oo behavior of the which gives the following expansion: 



7T.n T— >oo 
= 



gX ' 1 j + log(T/T ) + 



(4.3) 



The first term is simply, as already noted (Eq. ( |3.28| )) the free energy of a spin in the 
representation n x I, that is the impurity without any interaction with the electrons. This 
indicates a J = UV fixed point, with logarithmic corrections characteristic of asymptotic 
freedom. 

4-2. Low temperature (T < T , B <T). 

The region where both T and B are small compared to Tq corresponds to the vicinity 
of C = +oo. 

We shall in particular consider Fp for very low magnetic fields (B -C T); due to 
5 , C/(A^)-invariance the expansion takes the form (cf appendix A): 



1 N 



a=l 



which defines the magnetic susceptibility a] 1 . 

The extra term of Eq. ( |3.19| ) is now exponentially small, so that one obtains a sep- 
aration of the equations according to j < f,j = f,j > f. Note that to calculate the 
corrections to the ( = +oo value of x[\ we do need to take into account the shifts of ±%ix/N 
of (|3T9|) . 



• Underscreened case (/</). Again, it is easy to extend the ( — > +oo limit ( |3.24| ) to 
include a magnetic field; taking into account the power-law corrections, one finds that 



571 



loffY? t (e- B/T ) + 3 + 

10g ^-/V e ^ + log(T/T ^ 



(4.5) 



The first term is interpreted by saying that the impurity has been (under) screened by the 
electrons (cf fig. 2), and it now behaves like a spin in the n x (I — f) representation. This 



5 These power-law corrections are characteristic of a system of T.B.A. equations for which the 
index j takes an infinite set of values. 
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indicates a J = oo IR fixed point, the logarithmic corrections correspond to IR asymptotic 
freedom. 

In particular, since the electrons do not contribute to the zero temperature and zero 
magnetic field entropy (cf appendix B), we can extract the latter from Eq. Q4.5| ) and check 
that it is just the logarithm of the dimension of the representation of the ground state, as 
expected. 

For the susceptibility one finds 

nT ^o l n(N-n)(l-f)(l-f + N) 

1 ~ T N 2 ~l y ' ' 

which contains the normal low-temperature 1/T divergence. 

• Screened case (/ = /). The limit X/(C ~~ * +00) = 1 indicates that there is no term 
of order T in the expansion of FJ. This corresponds to the intuitive fact that the impurity 
has been completely screened by the electrons (J = 00 fixed point). Expanding around 
Q = +00 in the fusion equations for j = f one finds particularly simple (linear) equations 
and the expansion x 7 } = 1 + cs ^ sin(7m/iV)e _ ^ plus next corrections of the type e _fcl >, k 
integer; but not the constant in front of the first correction. Fortunately, for f = I, one 



can use an alternative method. One expands (|3.10| ) in the limit T <C Tq: 



f ttTo^ sin f J S 61 lfJ 

1 sin (ff ) 
To sin (f ) 



■T 



so that we are led to the computation of J 7 , which is done in appendix B. We only reproduce 
here the result (B.12): 



T ^o , mn(g) 
sin(f) 



n N 2 -1T 2 J_^(B a ) 2 



12 N + f T 4tt ^ T 



(4.1 



Furthermore, according to the form of the next corrections (e fc ^), has only integer 
powers of T in its expansion. 



The magnetic susceptibility cr^ is 



ff " T :°-l ^ ^ (4.9) 
f T sin(f) 2tt v ; 



1 sin(^) AT/ 
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which reaches a finite limit as T — > 0. This confirms that the impurity is completely 
screened in the IR. 

• Overscreened case (/ > I). This is the most interesting case. Independently of the 
magnetic field, one find the result ( |3.22| ), or more explicitly 



F? T ~°-T log 



' n N 

n n 

,p=l q=n-\-l 



sin(g — p + 1)9 
sm(q — p)9 



(4.10) 



where 9 = ir/(f + N). 

It is clear that here, the entropy at zero temperature S(T = 0) = log Xi(Gf) is n °t the 
logarithm of an integer number, a phenomenon which has a simple explanation (see next 
section). In particular, it is smaller than the naive value Xf-iO-SU(N)) one would obtain 
if one minimized only the interaction term in the Hamiltonian flli ]. This corresponds to 
the fact that the RG flow reaches a non-trivial IR fixed point J = J* characterized by 
Non-Fermi Liquid behavior. 

One should note that the number xf(Gf) involved is directly determined by the 
geometry of the fusion equations diagram (fig. 5). 



N-2 




f-2 



N-l O 



Fig. 5: Pictorial representation of the fusion equations j < f in the limit 
C — > oo, when they decouple from the equations j > f. 



The symmetries of the diagram imply equalities of zero temperature entropies for 
different impurities: obviously xT(Gf) = X^~ n (Gf) = Xf-i(Gf)- Furthermore, one can 
see on this graphical representation that the rank-level duality of the equations implies 
that the simultaneous interchange of n <-> I and N <-> / should leave xfi^f) unchanged. 
For example, this relates the cases of an overscreened impurity in a completely symmetric 
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representation (n = 1) and in a completely antisymmetric representation (1 = 1) (this 



particular duality can be found in the large iV limit in ||17|| ) 

Let us also calculate the next correction to the free energy. Since we are now dealing 
with a finite set of coupled equations (1 < j < f — 1, 1 < r < N — 1), one expects 
exponentially small corrections to the dominant behavior. We already know from ( [4.81 ) 
that x T f = 1 + 0(e~^); and * s connected to X r f-i by the fusion equations (even in 
the ( — > +oo limit: take ( |3.19| ) for j = f — 1), and from there on to all the Xai 3 < f ■ 
This implies that all the x T j have corrections of order e~^; but these may be subdominant 
compared to other exponentially small corrections. Hence, we try the ansatz 



X-(C)=X-(+oo) (l + aje-*) 



(4.11) 



for 1 < j < f — 1, with < r < 1. Expanding the fusion equations (|3.19|) around C — > +oo 
leads to the following set of linear equations for the a ^ : 



1 



« +1 +a 



1 + rjU+oo) 



+ 



1 



l + (7/r(+oo))- 



^(a r j+1 + a r j _ 1 ) = \a) 



where A = 2 cos(7rr/iV), and we have the boundary conditions a® = 



N 



(4.12) 



a'r 



a 



0. 



This is an eigenvalue problem; the dominant correction in ( [4.11] ) corresponds to the biggest 
A, so we are looking for the Perron-Frobenius eigenvalue and eigenvector!. The solution 



a r j = a sin(r<9) sin((iV - r)0) sin(j(9) sin((j + N)0) 
A = 2cos(2#) 



(4.13) 



satisfies ( |4.12|) , the boundary conditions and the positivity requirement > 0. Notice 
that ( |4.13| ) respects the rank-level duality of (|4.12| ) (exchange of r <-> j and N f). 

We finally find r = 2N/(f + N). We then need to treat separately the cases r < 1, 
r = 1, r > 1, which leads to the following discussion: 

★ For f<N,l<r<2 and the dominant correction is the e _< » correction; as the next 
correction coming from \ r f is °f order e _2< > , one can write 



-pn T— »0 rp 



(4.14) 



Note that the solution dj = 1 trivially satisfies ( 4.12 ) with A = 2, but not the boundary 
conditions. 
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where 7™ and a (the constant in af) are functions of the ratio B/T, so that 



a, oc 



N-f 



(4.15) 



The susceptibility tends to a finite value, as in the screened case, but the next correction 
has a non-integer power. 

* For / = N (the "self-dual", i.e. square diagram), r = 1: the exponential corrections 
of the x} are exactly of the same order as those of the Xj ? j ' < f, s ° that they have to 
be inserted as boundary conditions in (|4.12|) : a r j- now has a non-zero value which one can 
extract from (|4.8|) (this new boundary condition breaks the self-duality of the equations). 
In fact, there is no solution to this system of inhomogeneous linear equations, precisely 
because A = 2cos(7r/iV) (i.e. r = 1) is an eigenvalue of the matrix of this system. This is 
a "resonance" phenomenon, and as is usual in such degenerate cases we must replace the 
ansatz ( |4.11| ) with 

X -(C) = X, r (+oo) (1 + K + br jOe~ C ) (4.16) 

One can then find a solution to the corresponding inhomogeneous linear system with b 1 ^ of 
the type of ((CT3|): = b sm(Trr/N) sin(7rj /N) so that 



-T 



lo gx ?(G f ) + b??-log (jr) + 



(4.17) 



where again b is a S/T-dependent constant; 



n T^O 1 , (T 

ai " ^ l0§ [To 



(4.18) 



The susceptibility has a logarithmic divergence. 
* For / > N, t < 1, and we find: 



-T 



log X ?(G f ) + a? 



T 



+ . 



(4.19) 



(4.20) 



/-..\ 

The susceptibility is divergent as T — > 0, but with a non-integer power-law. 

In all three cases, the fact that the divergence of the susceptibility is always slower 
than 1/T indicates that the SU(N) spin of the impurity is completely screened at the IR 
fixed point (cf section 2). 
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5. Interpretation of the entropy at zero temperature. 

We shall first discuss in an abstract setting, how to calculate the entropic factor 
connected to a particle in a general theory. We shall then apply this to the impurity in our 
model, which will yield the zero temperature entropy (since the elctrons do not contribute 
to it). 

5.1. General principles. 

We consider a system in which states can be characterized by a quantum number R. 
As several states may have the same quantum number, we associate to each R a degeneracy 
cIr. Let us now imagine that we take a state which has the quantum number R, and add 
to it a particle: we obtain a new state which has the quantum number R' . Of course, 
for a given R and a given type of particle, not all R' are allowed; therefore we associate 
to the particle the adjacency matrix Arr> which is if the transition from R to R' is 
forbidden, and 1 (or in fact any positive integer in case of multiple possible transitions) if 
the transition is allowed. 

Let us now start from the vacuum, which conventionally has the quantum number 
R = 0. If we put one particle with associated adjacency matrix (Aq)rr/ into the system, 
the number of states allowed is by definition 



This expression is correct for a system of fixed finite length L; the entropy at zero temper- 
ature is equal to log O, that is the logarithm of an integer. However, we shall argue that 
this is in fact not the proper definition of the entropic factor associated to this particle in 
the thermodynamic limit. Indeed, as the size L of the system grows, at fixed temperature 
T, the average number of particles in the system necessarily diverges. In fact in a massless 
theory, it is obvious that the number of particles is of order TL ^> 1 (even if T is in the 
end sent to 0). 

Therefore, we come to the conclusion that to define the entropic factor of a given 
particle as L — > oo, one should always consider it as surrounded with a large number 
of other particles. One should then extract from the entropy of the resulting state the 
contribution of the particle we are interested in. 




(5.1) 



R 
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From this point of view the entropic factor is 

n = J2( A o A i 1A 2 2 ■ ■ ■ < p W* (5-2) 

R 

where we have introduced the numbers iVj and adjacency matrices A\ of the different types 
of particles. 

We shall make the assumption that all the matrices A{ commute. This will always be 
the case for us. Then in the limit iVj — > oo, we can write 

Ooc AoAfA^.-.A^ (5.3) 

where A^ is the biggest eigenvalue of Aj, which is necessarily real positive since A{ has 
positive entries, and which is associated with the Perron-Frobenius eigenvector common to 
all matrices Ai. The contribution to the entropy of our original particle has now become 
S = logA . 

An important remark is that the entropy does not depend on the degeneracy numbers 
d R ! This may sound slightly paradoxical, but it will in fact play a key role in what follows. 

Let us come back to our Kondo model. Here we have in fact two quantum numbers. 
We shall first introduce and study these two quantum numbers one at a time, constructing 
the adjacency matrices and computing their Perron-Frobenius eigenvalues according to the 
procedure described above. 

5.2. The SU(N) quantum number. 

Let us suppose first that the quantum number R is simply the SU (N) irreducible 
representation of the state. To a particle which belongs to the representation Rq we 
associate the adjacency matrix (A Rn ) RR ' which is the usual tensor product matrix in the 
space of representations; it is defined by the decomposition rule 

R ®R = Q(A Ro ) RR >R' 

R' 

i.e. matrix elements of A Ro are SU(N) Littlewood-Richardson coefficients. The degeneracy 
d R is the dimension of the representation R. Equation ( |5.1|) reads here 

n = J2( A Ro)<bRd R = d Ro (5.4) 

R 
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so that we immediately find the correct result that S = logO is the logarithm of the 
dimension of the representation Ro. Here, the formalism we developed is useless; it just 
gives in a complicated way the same result since the analogue of Eq. ( p.2| ) is 

n = ^RoKiKt ■ ■ ■ < p Md R = d Ro d N R \d N R i ...d N R i (5.5) 

R 

which immediately exhibits the Perron-Frobenius eigenvalues \ = d Ri of the adjacency 
matrices. 

5.3. The SU q (N) quantum number (restricted to level f). 

Let us first introduce the quantum group U q ($l(N)), with q = exp(2i6>) = exp(2z7r/ (/+ 
N)). We remind the reader that given a simple Lie algebra g, one can construct a family 
of deformations U q (g) of it: U q (g) is an algebra (in fact a Hopf algebra, see [|I8|] ) whose 
defining relations are deformations of those of g, with q the deformation parameter (when 
q — > 1 one recovers the undeformed universal enveloping algebra U(g)). Quantum groups 
are the natural symmetries of integrable models, even though the U q (sl(N) symmetry 
that we now consider is in a sense "hidden" in the Kondo model, and we rediscover it by 
computing thermodynamic quantities. 

Let us turn to the representation theory of U q (sl(N)). For generic q, it is identical 
to that of sl(N) (and irreducible representations can be depicted e.g. by Young tableaux 
as in appendix A). However, for q a root of unityH, the situation is more complicated: 
a subclass of irreducible representations, the "good" representations behave like for any 
value of q; but the others have a more complicated behavior (in particular several Young 
tableaux merge together into indecomposable but not irreducible representations). In our 
case, when q is a primitive (/ + iV) th root of unity, the "good" representations, which we 
also call u SU q (N) representations" for simplicity, are characterized by the property that 
the number of columns of their Young tableau n 1 + ■ ■ ■ + n N ~ 1 = m 1 —m N < f. We denote 
their set by P+(N,f)E. 

Fortunately, it can be shown that one can consistently restrict the representation 
theory U q (sl(N)) to the "good" representations. This is the situation we consider now: 
the SU q (N) quantum number is precisely an element of the set P + (N : f). The justification 

7 Here, we consider the so-called "restricted specialization" of U q (sl(N)) 

8 We borrow this notation from affine Lie algebra theory, since P + (N,f) can also be defined 
as the set of integral dominant weights of ^4jv-i a ^ level /. 
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of this choice is that, from the form of the B.A.E., we conjecture that in the spin sector of 
the Kondo model, only such representations are present. In fact, this can also be guessed 
from the similar situation in the WZW SU(N) level / theory. 

Let us comment on the physical implication of such a quantum number. One-particle 
states are now organized in multiplets of SU q (N) with certain "good" representations, 
(for us, they are the fundamental representations of SU q (N)); however, as soon as one 
considers multi-particle states, it becomes clear that the Hilbert space of our theory does 
not have a Fock structure any more: indeed by performing the tensor product of enough 
representations of P + (N : f) one necessarily obtains representations with more than / 
columns and which therefore do no belong to P + (N, f). A truncation is needed to keep 
only the "good" representations in the tensor product. 

The non-Fock structure of the states is most easily imagined by considering the phys- 
ical particles as kinks ||19[| , that is states interpolating between different classical vacua 
labeled by the SU q (N) quantum number. This is the same "kink structure" which ap- 



pears, for example, in RSOS models ||20|| . In the end, this non-Fock structure will be 
responsible for the fractional zero temperature entropy. 

This truncation (and the associated kink structure) is apparent when one builds the 
adjacency matrices (A t a)rri\ they are the truncated tensor product matrices for the quan- 
tum group U q (sl(N)) associated to the rectangular Young tableaux r x j. In the kink 
picture, the adjacency matrices describe the possible transitions between the different 
classical vacua. Here we do not give an abstract algebraic construction for A r - but rather 
a simple explicit procedure. 

We start with A\ which is obtained by taking the corresponding standard (i.e. untrun- 
cated) tensor product matrix and restricting it (in the naive sense) to P+(N, f). Explicitly, 
this means that (A\)rr/ = 1 if there is a a G {1, . . . , N} such that R' = R+e a , otherwise. 
We have shown in fig. 6 the set -P+(3, 3); each arrow corresponds to the addition of a e a , 
a = 1 ... 3, so that A\ is precisely the adjacency matrix of the resulting graph. 

This way we have defined the matrix A\. We could use the same restriction procedure 
for any matrix A\ corresponding to a fundamental representation of SU(N); explicitly, 
{A[)rr> = 1 if there are ai, . . ., a r all distinct such that R' = R + e ai + ■ ■ ■ + e Qr , 
otherwise. Fig. 7 gives an example of such a construction. 
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Fig. 6: The set -P+(3, 3) and the graph of A\. 




□ m 

Fig. 7: The set i~+(4, 2) and the graph of A\, which consists of the edges 
with double arrows. 

The higher matrices A r -, j > 1 cannot be obtained in the same way; one must use 
the fusion procedure. In the case of rectangular Young tableaux, the matrices A r - can be 
defined by imposing the by now familiar fusion equations 

A^A) = A] +1 A]~ l + A r j+l A)_ x (5.6) 

with the boundary conditions A® = = Aq = 1. It should be noted that the fusion pro- 
cedure requires only A\ to build all the other Ay, however, the sole use of eqs. ( |5 . 6| ) requires 
all the A\, so as to have proper initial conditions to the recursion A r - = f(A < j_ 11 A q -_ 2 ). In 
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fact one has the following explicit solution in determinant form: 

A r j = det(Ar a+6 )i< , 6 < i (5.7) 

(where A\ = for r < or r > N). 

One can finally check that the largest eigenvalue of Aj, which like all eigenvalues is a 
SU(N) character, is precisely x^(G/)i. 

5.4- The Kondo model. 

In the Kondo model, as already mentioned, the two quantum numbers described above 
naturally appear in the Bethe Ansatz description through the appearance of the ^-strings 
for 1 < j < / and j > f. We shall now try to apply the general method outlined above to 
the impurity in interaction with the electrons, i.e. the (under, over) screened impurity. By 
inspection of the ground state (section 2), we can infer that the impurity never has both 
adjacency matrices non-trivial: for the underscreened case only the unrestricted SU(N) 
adjacency matrix is non-trivial, while the SU q (N) adjacency matrix is the identity matrix, 
and the converse statement for the overscreened case. This means that we can hide one 
of the two quantum numbers inside the degeneracies cLr of the other; as we have already 
noted that these numbers d,R play no role in the thermodynamic limit, we can blindly 
apply the results of the two previous subsections: 

• Underscreened case (/ < /). We apply the result of section 5.2 to the impurity, 
which we assume belongs to the representation n x (I — /) due to the screening of the 
electrons. Our procedure gives the correct result that the entropy at zero temperature is 
the logarithm of the dimension of this representation. 

• Overscreened case (/ > /). We assume that the adjacency matrix corresponding 
to the impurity is the SU q (N) truncated tensor product matrix A^~ n ; this non-trivial 
adjacency matrix is due to the (N — n) x I electrons (or more precisely their spin sector 
component) on the site of the impurity, which are screening it: indeed to completely screen 
the spin of the impurity (as was found in section 2) one needs precisely (N—n) x / electrons. 

In this case too, one finds the correct result that the entropy is the logarithm of 
xr n (Gf) = X?(G f ). 

9 In the quantum group terminology this number is called the quantum dimension of the 
representation r x j. 
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Appendix A. SU(N) representations and characters 

An irreducible representation (irrep) R of U(N) is characterized by its highest weights, 
which form a sequence of N integers m 1 > m 2 > . . . > m . m a is the length of the a th 
row of the corresponding Young tableau (fig. 8). 
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Fig. 8: A generic Young tableau. 
We shall write: 

N 

R = ^m a e a (A.l) 

a=l 

As we are dealing with SU(N) (and not U(N)) irreps, we simply add the rule that two 
irreps whose highest weights only differ by a constant are equivalent. This amounts to 
imposing the relation: e 1 + e 2 + • • • + e N = 0. One possible convention is then to suppose 
that m N = 0, the other highest weights m a being non-negative (as on figs. 6 and 7). 

One can also introduce another set of numbers to describe R: one defines n r to be the 
number of columns of given size r, 1 < r < N — 1. This corresponds to the decomposition 

N-l 

R=^n r uj r (A.2) 

r=l 

where u r = X/I=i e<1 ^ s ^ ne rth fundamental weight of SU(N). We have n r = m r — m r+1 . 
In the case of a Bethe Ansatz state, using Eq. ( |1.21|) , one can show that 

n r = n r -2j2C qr M q (A.3) 
q 
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where n r is the contribution of #0, and 2C qr = 2C qr (n = 0) = 25 qr - (S qr+1 + S^' 1 ) is 
the usual Cartan matrix of An-i- Rewritten explicitly, this gives Eq. (|2.6| ). 

To a representation R and an element G of SU(N) (or, by analytic continuation, of 
SL(N)) we associate the character Xr(G) defined as the trace of G in the representation 
R. It satisfies the basic SU (N) invariance property: 

Xr(G) =x fl (OGO- 1 ) \/neSU(N) 

which implies that it only depends on the eigenvalues z a (1 < a < N) of G. For irreps, 
one has the following explicit formula: 

det ^ z a^m b +N-b^ 



= ^ml^f (A ' 4) 



where det((z a ) N ~ b ) = A is simply the Van der Monde determinant of the z a . 

Let us now suppose that the z a are real positive, so that one can order them: z 1 > 
z 2 > . . . > z N . Let us furthermore assume that the representation R becomes large, in 
the sense that the n r 3> 1. Then one can replace ( |A.4j) with the following estimate: 



1 N 

Xr(G) ~ ^H(z a r a+N ~ a (A.5) 
0=1 

Setting z a = exp(—B a /T) and expressing the m a in terms of the M r , one obtains Eq. 

For a large rectangular Young tableau r X j, we have n 9 = j'5 9r so that only n r ^> 1. 
A more careful estimate gives 

x r {G )^^f[(z«r a+N - a (A.6) 

a=l 

where Ai (resp. A2) is the Van der Monde of the z a for 1 < a < r (resp. r + 1 < a < N). 
It is quite close to ( |A.5|) in the sense that the character is still approximated by a constant 
times the highest weight vector contribution. 

More explicitly, for z a = exp(—B a /T), we have 

rr -B/T\i^°° -jT'' B a /T /» „n 

Xj\ e ' ) oc e J ^a=i 1 (A. 7) 



(compare with Eq. ( |3.32| )) . 
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Let us finally expand Xii( e ~ B ^ T ) f° r -B <C T, i.e. around XrO-su(N))- Due to the 
S'£7(iV)-invariance, the linear term is necessarily proportional to the only linear invari- 
ant ^2^ = iB a , which in our case is zero. In the same way, there are only two quadratic 
invariants, (J2a=i S ") 2 = and E^LiO 8 ") 2 , so that 

N 

\og X R(e- B/T ) = log XR (lsu(N)) + -7T E( Sa /^) 2 + • • • (A.8) 

a=l 

This motivates the definition of the magnetic susceptibility given in (JOl) . To compute the 
constant C, one uses the well-known identity: A\ r — C2(R)xr where A is the Laplacian on 
SU(N) and C^-R) is the quadratic Casimir of the irrep R. With the usual normalizations 
this means that C = C 2 (R)/(N 2 - 1). 



Appendix B. Computation of the free energy of the electrons 

This appendix presents the computation (from the Bethe Ansatz), at given tempera- 
ture T and magnetic field B, of the free energy per unit length of the electrons 



T= J rfAG rl *s(A-l)log(l+^(A)) 



(B.l) 



We shall take the scaling limit only at the end of the calculations. 

The key observation is that 2nD5jjG rl *s(A-l) = (d/dA)^ (cf replacing 
g r - with the l.h.s. of (|3.7|) , and inserting this into (|B.1|) gives 



^ = -^E/ dK JR (log(l + (^)- 1 )-E G9r *^^ 1 °g(l + ^) )log(l + i£(A)) 

r,j \ q,k J 

(B.2) 

The double integral which appears in ( B.2|) 

/= J2 /^AdA'^(A-AOlog(l + ^(A))log(l + ^(A')) (B.3) 

j,k,q,r 

where A q - r k = (d/dA)(G qr *Cjk), would seem to vanish due to the change of sign under the 
simultaneous transformations q r, j <-> A;, A «-» —A. In fact it does not, because of the 
non-uniform convergence of this integral as A, A' — > ±00, and the correct result is: 

+00 



I= \ Y, Gi r (K = 0)C jk (K = 0) 



j,k,q,r 



log(l + rjj) log(l + vl) 



(B.4) 
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We now use Eq. ( |3.1| ) in the limit A — * ±00 

n r (±Oo) 

log(l + »£(±oo)) - Cqr ^ = °) G j^ = 0) log(l + (^(ioo))" 1 ) = ^— -L (B.5) 



q,k 



to express in a more suggestive form /: 



log(l + rfl) log(l + (77 



+00 



+ 5 £ G« r (« = 0)C ifc (/c = 0) 



3,k,q,r 

The free energy then takes the form 
1 



f log(l + ^) 



+00 



(B.6) 



T 2 
9tt 



[d(log(l + (tjJ)" 1 )) log(l + 17J) + d(log(l + tjJ)) log(l + (tjJ)" 1 )] 



j,k,q,r 



g r Aog(l + V l) 



+00 



(B.7) 

(k = is implied for all kernels now). We shall compute this expression in the scaling 
limit. 

The first term is usually written in terms of Rogers' dilogarithmic function L, by 
taking 1/(1 + 77) as new variable; in fact one can easily show that, due to cancellations 
between A = +00 and A = —00, the quantity between the brackets is precisely equal to 



the following (known) dilogarithm sum EI 



/ N-l 



4^ ^ I 1 + 77^1 



3=1 r=l 



+ r]U+oo) J 6 N + f 



(N 2 - 1) 



(B.8) 



Concerning the second term, one might be tempted to replace ^ r ■ ; G qr Cjk <7j(±oo) 
with its value (^(±00) and then evaluate the sum, which turns out to be zero. This is in 
fact incorrect because of the divergences which occur at j = 00 due to the presence of the 
magnetic field B. The correct procedure is to isolate the term where B appears 



+00 j 



'og(l + C 



(B.9) 



q,r 



k=l 



3=1 
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where we have introduced a cutoff J for the last sum (which would naively vanish as 
J — > oo). Using Y, J j=ij C jk = 8 k j(J + l)/2- 5 kJ+1 J/2 and (|3.13|) we find: 



rp r 1 +°° 



T B = — lim V6 r 



(J + l)log X }- Jlogx} +1 



(B.10) 



The explicit expressions oo) = Xj( e S/,T ) 5 = Xj-/( e B//T ) (j > /), and the 

estimate ( A..7) allow to compute ( |B.1U| ); the result is: 



^ = -^5> a ) 2 ( B - n ) 

a=l 

Putting (|B.8| ) and ( [B.llp together, we obtain the final expression for the free energy: 



F=-T 271 f 



1 - 

( N2 - 1 )-JzfY,( Ba ) 2 (B.12) 



12 N + f ' 4tt 

J a=l 

Let us comment on this result. The first term is the central charge term for a theory 
with only one chirality and central charge c = (iV 2 — l)f/(N + /): this is what we expect 
since the spin sector of the electrons is a WZW SU(N) level / theory. Of course once one 
puts together the spin, flavor and charge sectors, one recovers the total central charge 

c=(N 2 - + (f ~ Vj^Tf +l=Nf (B.13) 

for iV / free complex fermions. 

The second term can be directly obtained as the free energy of the electrons in the 
magnetic field B: indeed the flavor and charge sectors do not contribute to it. One has 
the simple formula 

Tb = ~Tf f °° f- J2 log(l + e-^/ T ) (B.14) 

J - U a=l 

where T> is some momentum cutoff of the order of D. The factor of / comes from the / 
flavors. After sending T> to infinity and discarding the divergent T-independent part, one 
obtains flBTnp . 
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